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Measurements of D and D s decay constants at CLEO 
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Using CLEO data collected at 3370 MeV and 4170 MeV, we determine f D+ = (205.8 ± 8.5 ± 2.5) MeV and an interim 
preliminary value of f D + = (267.9 ± 8.2 ± 3.9) MeV, where both results are radiatively corrected. They agree with the 
recent most precise unquenched Lattice-QCD calculation for the D + , but do not for the . Several consequences 
are discussed. 

1. Introduction 

Leptonic decay — » l^v is described by the annihilation of the initial quark-antiquark pair into a virtual W + 
that materializes as a pair (Fig. 1). 




Figure 1: Decay diagram for Df, — > l + v. 



The decay rate is give by 

r^J) - ^) - ^f Dt y e m Dt) (i - JjL-j \v cd(s) \\ (1) 

where f D + is decay constant [1], related to the overlap of the heavy and light quark wave- function at zero spatial 
separation, Gf is the Fermi constant, m D + is the mass, me is the final state charged-lepton mass, and V^( s ) is 
a CKM matrix element, taken as V c d — V us = 0.2256 and V cs = V u d = 0.9742. Thus, measurement of purely leptonic 
decays allow a determination of the decay constant f n + . 

00 

Meson decay constants in the B system are used to translate measurements of BB mixing to CKM matrix elements. 
Currently, it is not possible to determine fg accurately from leptonic B decays, so theoretical calculations of fg are 
used. Since the B® meson does not have l + v decays, it will never be possible to determine fg s experimentally, thus 
theory must be relied on. If calculations disagree on D mesons, they may be questionable on B mesons. If, on the 
other hand new physics is present, we need understand how it effects SM based predictions of the B decay constants. 
Decay constants can be calculated using lattice quantum-chromodynamics (LQCD). Recent calculation from Follana 
et al. using an unquenched LQCD predicts f D + = (207 ± 4) MeV and f D + = (241 ± 3) MeV [2]. 

We use the reactions e+e~ -> D~D+, and e+e~ -> D*~D+ or D-D*+. The D+ is studied at 3770 MeV using 
818 pb" 1 of data [3]. And the D+ is studied at 4170 MeV, using 424 pb" 1 for (i+v and t+v, t+ -> tt+P, and 300 
pb _1 for t + v, r + — > e + vv. (Eventually CLEO will present results using 600 pb _1 .) 



2. D+ -> £+ 



v 



We fully reconstruct D [4] as a tag and examine the properties of the other D + , which can be found even if 
there is a missing neutrino in the final state. This method is called as "double tag" technique. To reconstruct D~ 
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tags we require that the tag candidates have a measured energy consistent with the beam energy, and have a "beam 
constrained mass", tobc, consistent with the D~ nominal mass [5], where mec = \J E^ cam — Pi) 2 , E\> cam is the 
beam energy and i runs over all the final state particles three-momenta. Fig. 2 shows the uibc distribution summed 
over all the decay modes we use for tagging. Selecting events in the mass peak we count 460,055±787 signal events 
over a background of 89,472 events. 
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Figure 2: The beam-constrained mass distributions summed over D~ decay candidates in the final states: K + -k~it~, 
K + tt~iy~tt°, Ksit~ , Kstt~iy~tt + , Kstt~iv and K + K~1T~. 

We then search for signal events with one and only one additional track with opposite sign of charge to the tag, 
not identified as kaon. The track must make an angle > 25.8° with respect to the beam line (90% of the solid angle), 
and in addition we require that there not be any photon detected in the calorimeter with energy greater than 250 
MeV. The latter selection can highly suppress D + — > 7r + 7r° background. We separate these events into two cases, 
where case (i) refers to muon candidates depositing < 300 MeV, characteristic of 98.8% of muons and case (ii) is for 
candidates depositing > 300 MeV, characteristic of 45% of the pions, those that happen to interact in the calorimeter 
and deposit significant energy. 

We look for D + — > [i + v by computing the square of the missing mass 

MM 2 = (£ bcam -E^) 2 - (-p D - - p M+ ) 2 , (2) 

where Pd- is the three-momentum of the fully reconstructed D~ , and i? M +(p M +) is the energy (momentum) of the 
[i + candidate. The signal peaks at zero for n + v and is smeared toward more positive values for t + v, t+ — ► tt + ~F. 

The fit to the case (i) MM 2 distribution shown in Fig. 3 contains separate shapes for signal, 7r + 7r°, K ir + , t + v 
(t + — > 7r + P), and a background shape describing three-body decays. Here we assume the SM ratio of 2.65 for the 
ratio of the t + v/^ + v component and constrain the area ratio of these components to the product of 2.65 with 
B(t+ -> tt+P) = (10.90±0.07)% [5] and the 55% probability that the pion deposits <300 MeV in the calorimeter. The 
7T + 7r° background are fixed at 9.2 events obtained from Monte Carlo (MC) simulation. The normalizations of the 
signal, K tt + , and 3-body background are allowed to float. The K ir + shape is obtained from double tag events of 
D — > K ~7r + , D° — > K + ir~ where we ignore one kaon to calculate the MM 2 . All other shapes are obtained from the 
MC simulation. 

The fit yields 149.7±12.0 /j, + v signal events and 25.8 t + v, t + — > tt + P events (for the entire MM 2 range). We also 
perform the fit allowing the t + v, t + — > Tt + i> component to float. Then we find 153.9±13.5 [i^v events and 13.5±15.3 
t + v, t + — > tt + v events, compared with the 25.8 we expect in the SM. Performing the fit in this manner gives a result 
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Figure 3: Fit to the MM 2 for case (i). The points with error bars show the data. The black (dashed) curve centered at zero 
shows the signal fi + u events. The dot-dashed (red) curve that peaks around 0.05 GeV 2 shows the D + -> r + 
component. The solid (blue) Gaussian shaped curve centered on the pion-mass squared shows the residual 7r + 7r° component. 
The dashed (purple) curve that falls to zero around 0.03 GeV 2 is the sum of all the other background components, except the 
tail which is shown by the long-dashed (green) curve that peaks up at 0.25 GeV 2 . The solid (black) curve is the sum 
of all the components. 

that is independent of the SM expectation of the D + — > t + v rate. To extract a branching fraction, in either case, we 
subtract off 2.4±1.0 events found from simulations and other studies to be additional backgrounds, not taken into 
account by the fit. 

We find B(D+ -> n+u) = (3.82 ± 0.32 ± 0.09) x 10~ 4 . The decay constant f D + is then obtained from Eq. (1) using 
1040±7 fs as the D + lifetime [5] and 0.2256 as \V c d\. Our final result, including radiative corrections is 

f D + = (205.8 ± 8.5 ± 2.5) MeV . (3) 

3. d+ -> e+u 

s 

To reconstruct the tag, the difference here than the D + case is that we need to include an additional photon from 
D* in the tag. The D~ tags we reconstructed can either from directly produced Dj mesons or those that result from 
the decay of D* mesons. We calculate the missing mass squared MM* 2 recoiling against the photon and the D~ tag, 

MM* 2 = {Ecu -E Ds -E 7 ) 2 - (pcm - PD S - P 7 ) 2 , (4) 

here -Ecm (Pcm) is the center-of-mass energy (momentum), Ejj s (jPd 3 ) is the energy (momentum) of the fully 
reconstructed D~ tag, and E 7 (p 7 ) is the energy (momentum) of the additional photon. We determine number 
of tags by simultaneously fit to the invariant mass (Md s ) and MM* 2 , shown in Fig. 3. The signal is fit to a sum 
of two Gaussians for the Md 3 and a Crystal Ball function for the MM* 2 . The tail parameters of the Crystal Ball 
function are obtained from fully reconstructed D*D S events. The background has two components: either comes 
from the background under the invariant mass peak (fake D s ), or is due to random photon combinations. The former 
background is modeled by linear function for the mass and 5th order Chebyshev polynomial for the MM* 2 . Since the 
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latter background is from a true D s , its mass distribution has the same shape as that from the signal, while its MM* 2 
is modeled by another 5th order Chebyshev polynomial function. The total number of single tags is 30848 ± 695 ± 925 
in the invariant mass signal region (±17.5 MeV from the nominal D s mass) and MM* 2 e[3.782, 4.0] GeV 2 . 




Figure 4: Distributions of the invariant mass (left) in MM* 2 G [3.5, 4.25] GeV 2 and MM* 2 (right) in the invariant mass 
\Md s — 1968.3 MeV| < 17.5 MeV in the final states: K + K~ix~ , KsK~, r]-K~; -q — ► 77, r)'n~; 7/ — > ir + ir~r], r\ — > 77, <j)p~; 
<p — > K + K~ , p~ — > 7T~7r°, 7r + 7r~7r~, K*~K*°\ K*~ — > Kgiv~ , K*° —t K + tv~, t]p~ \ 77 — > 77, p~ — > 7r~7r°, and r/'ir~; 
rj —* 7r + 7r~7. Data (point) and fit function (red curve) are shown. The dashed (blue) curve corresponds to the fake D s 
background, and the dashed (green) curve to the random photon background. 



Similarly as the D + , we reconstruct the signal side. We veto events with an extra neutral energy cluster > 300 
MeV (it is 250 MeV in the D + case). It is highly effective in reducing backgrounds, especially for Df — > ir + ir°, rjir + 
and the processes D^*>D^*'. The missing mass squared, MM 2 , evaluated by taking into account the observed 
Dj , and photon should peak at zero; 

MM 2 = (E CM - E Ds -E y - E^f - (p CM - Pd 3 - P 7 - P M ) 2 - (5) 

We also make use of a set of kinematical constraints and fit each event to two hypotheses one of which is that the 
tag is the daughter of a D*~ and the other that the D* + decays into with the Df subsequently decaying 

into [i + v. In addition, we constrain the invariant mass of the Dj tag to the know D s mass. This gives us a total of 7 
constraints. The missing neutrino four-vector needs to be determined, so we are left with a three-constraint fit. We 
perform a standard iterative fit minimizing x 2 . We choose the fitted MM 2 from the hypothesis giving the smaller 
X 2 . The MM 2 distribution from data is show in Fig. 5. After fixing the ratio of t + v/ n + v to the SM value we find 
f D + = (268.2 ± 9.6 ± 4.4) MeV. 

We can also use the decay mode r + — > e + vv. This result has already been published. [6] The technique here 
is to use only three tagging modes: 07r~, K~K*° and KgK~ , to ensure that the tags are extremely clean. Then 
events with an identified e + and no other charged tracks are selected. Any energy not associated with the tag decay 
products is tabulated. Those events with small extra energy below 400 MeV are mostly pure Df — » t + v events. 
After correcting for efficiencies and residual backgrounds we find f D + = (273 ± 16 ± 8) MeV. 



4. Conclusions 



The preliminary CLEO average is f D + = (267.9 ± 8.2 ± 3.9) MeV (radiatively corrected). Averaging in the Belle 
result [7] f D + — (269.6 ± 8.3) MeV, which differs from the Follana et al. calculation [2] by 3.2 standard deviations, 
while the result for f D + — (205.8 ± 8.5 ± 2.5) MeV is in good agreement. This discrepancy could be due to physics 
beyond the standard model [8], or systematic uncertainties that are not understood in the LQCD calculation or 
the experimental measurements, or unlikely statistical fluctuations in the experimental measurements or the LQCD 
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Figure 5: The MM 2 distribution. The dashed (grey) Gaussian shaped curve peaked near zero is the fi + u component, while 
the dashed (purple) curve that rises sharply from zero and then flattens out shows the t + v component. The two lines are 
background components. The solid curve shows the sum. 



calculation. Fits to the CKM matrix parameters use theoretical predictions of Jb s / B d ■ As similar calculations are 
used for fsj B d , we need to be concerned with them. 
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